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DYNAMIC  TRANSVERSE  LOADING  OP  PEAMS 
OP  A MATERIAL  EXHIBITING  LINEAR  STRAIN-HARDENING1 

by 

Margaret  F.  Conroy2 

Abstract.  Tho  object  of  this  report  is  to  show  that  in  certain 
circumstances  tho  plastic  deformation  of  beams^  made  of  a material 
with  linear  strain-hardening  and  subject  to  dynamic  transverse 
loading,  can  be  determined  by  the  techniques  usod  in  solving 
elastic  problems.  In  particular,  the  differential  equation  of 
motion  for  such  beams  is  in  some  instances  of  the  same  form  as 
in  the  corresponding  linear  elastic  case,  and  so  any  of  tho 
methods  employed  for  solving  elastic  beam  problems,  such  as  the 
normal  mode  method,  Ieplace  transform  mothod,  or  Boussinosq!s 
solutions  for  infinite  beams,  can  bo  used.  Because  of  this 
linear  character  of  the  differential  equation  of  motion  en- 
countered in  the  analysis  presented  here,  it  is  also  shown  that 
some  initial  motion  problems  for  beams  undergoing  large  plastic 
deformations  due  to  transverse  loading  can  bo  solved  by  supor- 
posing  solutions.  In  these  problems  tho  disturbance  part  of  the 

solution  is  obtained  by  somo  elasticity  technique  and  is  then 

\ 

supcrposcc  on  the  initial  motion  of  the  beam. 

The  method  of  solution  is  demonstrated  by  moans  of  sovoral 
examples  involving  finite  beams.  The  first  oxample  is  an  initial 
motion  problem  and  illustrates  tho  mothod  of  superposition.  The 

1.  Tho  results  presented  in  this  paper  were  obtainod  in  the  course 
of  research  sponsored  by  the  Office  of  Naval  Research  undor 
Contract  N7onr-35801  with  Brown  University. 

2.  Researoh  Assistant,  Graduato  Division  of  Appliod  Mathematics, 
Brown  University. 
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disturbance  part  of  tho  solution  is  found  by  the  normal  modo 
technique,  commonly  used  to  solve  forcod  vibration  problems  in 
elasticity*  An  approximation  to  the  solution  of  tho  same  problem 
is  found  by  means  of  Boussinosq»s  solution  for  an  infinite 
clastic  beam.  Tho  socond  oxamplo  is  an  initial  stress  problom 
for  a simply  supported  beam.  Again  tho  method  of  superposition 
is  used*  The  last  oxamplo  is  a free  boundary  problom  for  a 
cantilever  beam.  This  problom  is  solved  by  an  inverse  method 
whereby  the  form  of  tho  solution  is  assumed  and  tho  physical 
problom  associated  with  this  solution  is  then  determined*  (A 
similar  typo  of  analysis  is  used  in  finding  solutions  of  elastic 
plane  strain  problems  by  the  consideration  of  simple  polynomial 
solutions  of  the  biharmonic  equation*) 

The  general  problem  of  the  determination  of  tho  plastio 
deformation  of  beams  subject  to  dynamic  transverse  loading  is 
very  difficult  to  handlo*  Thus  far  only  ono  such  problem  has 
been  solved  for  a beam  mado  of  strain-hardening  material*  H.  F. 
Bohnonblust  [l]*  has  made  an  elastic-plastic  analysis  of  tho 
problem  of  an  infinite  beam  subject  to  a constant  velocity  im- 
pact. In  view,  then,  of  tho  difficulty  encountered  in  solving 
problems  of  this  type,  tho  analysis  carried  out  in  this  report, 
while  it  applies  only  to  a very  special  class  of  probloms,  sooms 
woll  worthwhile,  since  it  docs  add  a group  of  tractable  solu- 
tions to  the  litcraturo.  It  is  believed  that  this  analysis 
represents  the  first  treatment  of  the  plastic  deformation  of 
finite  beams,  made  of  a strain-hardening  material  and  subjoct 
to  dynamic  transverse  loading* 

* Numbers  in  square  brackets  refer  to  the  Bibliography  at  the  end 
of  the  paper* 
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x Distance  along  the  beam, 
y Deflection. 

K fi  y/fix  , the  curvature  of  the  deflection  curve, 
m Mass  per  unit  length  of  the  beam. 

M Bending  moment. 

Q Shearing  force, 

t Time. 

M0  The  limiting  value  of  M for  rigid  body  motion. 
Foot-pound-second  units  are  used  throughout  this  report. 

1.  Introduction  and  basic  assumptions.  The  analysis  carried  out 
in  this  report  is  introduced  to  determine  the  permanent  deforma- 
tion of  beams  subjected  to  transverse  loading  of  such  a magnitude 
that  the  plastic  strains  produced  are  large  compared  to  tho 
elastic  strains.  It  is  shown  here  that,  in  some  instances,  a 
satisfactory  approximation  to  t he  solution  of  such  problems  oan 
be  obtained  by  means  of  the  same  tochniques  used  in  elastio  beam 
theory. 

The  analysis  is  based  on  tho  assumption  of  a linear  strain- 
hardening  bending  moment -curvaturo  relationship  of  the  typo  shown 
in  Fig.  1.  Thus,  o la stic  strains  are  ignored  and  the  boam  is 
considered  to  be  either  rigid  or  plastio.  An  approximation  for 
the  solution  of  tho  actual  olastic-plastic  problem  is,  then, 
obtained  by  neglecting  the  elastic  strains  and  carrying  out  a 
plastic-rigid  typo  of  analysis.  It  is  expected  that  this  approxi- 
mation will  be  satisfactory  when  the  plastic  strains  are  largo 
compared  to  tho  elastic  strains. 
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Fig.  1 


The  analysis  is  also  based  on  the  assumption  that  the  rate 
of  change  of  curvature  of  the  beam  is  everywhere  greater  than 
or  equal  to  zero,  or  everywhere  less  than  or  equal  to  zero,  for 
the  range  of  time  considered,  and  this  of  course  places  certain 
restrictions  on  the  initial  conditions,  and  boundary  conditions 
for  application  of  this  analysis,  as  described  below, 

Undor  these  assumptions  the  bending  moment-curvature  rela- 
tionship of  the  beam  for  plastic  flow  is  the  linear  relationship 

2 

M = + M0  + bK  = + M0  + b (1) 

where  the  sign  of  MQ  is  the  same  as  the  sign  of  the  rate  of 
change  of  curvature  of  the  team.  The  differential  equation  of 
motion  for  the  beam,  which  is  obtained  from  tho  equilibrium 
conditions 

Q = ®H  , + m = 0 

ax  ’ dx  712 


> 
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becomes 


Ox*  c2  8t2 


(2) 


where  -L  * -SI  . 
c2  b 

This  equation  of  motion  is  of  the  same  form  as  in  the  elas- 
tic beam  theory,  where  the  flexural  rigidity  has  been  replaced 
by  b.  Thus,  under  the  assumptions  stated,  it  is  clear  that  the 
plastic  deformation  of  a beam  can  be  determined  by  the  same 
methods  used  to  solve  the  corresponding  elastic  problem  and, 
since  the  equation  of  motion  is  linear,  superposition  of  effects 
is  permissible. 

The  assumptions  stated  here  restrict  the  application  of  the 
analysis  to  a very  special  class  of  problems*  However,  as  men- 
tioned before,  the  general  problem  of  the  plastic  deformation  of 
beams  subject  to  dynamic  transverse  loading  has  not  yielded  to 
theoretical  analysis,  hence  the  analysis  given  here  represents 
definite  progress* 


2.  Problem  I — An  Initial  motion  problem*  Fcr  the  application 
of  this  method  to  an  initial  motion  problem,  the  initial  motion 
of  the  beam  must  be  such  that  the  final  rate  of  change  of  curva- 
ture of  the  beam  after  loading  is  everywhere  greater  than  or 
equal  to  zero  or  everywhere  loss  than  or  equal  to  zero*  That  is, 
the  rate  of  change  of  curvature  of  the  beam  due  to  the  initial 
motion  must  be  such  as  to  eliminate  any  oscillations  in  the  sign 
of  the  rate  of  change  of  curvature  due  to  the  transverse  loading. 
An  example  of  such  an  initial  motion  is  afforded  by  the  following 
problem* 
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Let  a cantilever  beam  of  uniform  cross  section  and  longth  -t 
have  ono  end  (x  = 0)  fixed,  and  bo  subject  to  an  initial  motion 

y * 2.2?tac2  - si  (3) 

which  is  maintained  for  all  time  by  a bonding  momont,  MQ  + b(4«5t 
“ •jjjj ) , applied  at  tho  free  end  of  the  beam.  For  t > 0 a bonding 
moment  9bt  is  appliod  to  tho  froo  end  (x  - t)  of  tho  boamc  Tho 
problem  is  to  dotormine  tho  subsequont  motion  of  tho  beam* 

The  initial  motion  and  tho  appliod  bonding  moment  given  hero 
woro  choson  rather  arbitrarily.  Any  othor  choices  would  servo 
as  well,  providod  thoy  aro  such  that  tho  rato  of  chango  of  curva- 
ture of  the  beam  after  loading  is  ovorywhero  greater  than  or 
equAl  to  zero  (or  everywhere  loss  than  or  equal  to  zoro). 

Due  to  tho  prescribed  Initial  motion,  tho  initial  rate 
of  -ohange  of  curvature  of  the  boam  is  evorywhoro  greater  than 
zero.  Thus,  the  bonding  momont-curvaturo  relationship  is 

Mj  = Mq  + bKj  (4) 

whore  Kj  and  Mj  are  tho  bonding  moment  and  curvature  correspond- 
ing to  tho  initial  motion.  If  tho  final  rato  of  chango  of  curva- 
ture of  the  beam  duo  to  tho  initial  motion  plus  the  disturbance 
remains  ovorywhero  positive,  tho  linear  bending  moment-curvature 
relationship  of  Eq.  (1)  remains  valid  and  from  Eqs.  (1)  and  ( f) 
it  is  clear  that 

Md  = «D 

where  KD  and  MD  aro  the  additional  bonding  moment  and  curvature 
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The  final  deformation  of  the  beam  can  be  obtained  by  super- 
posing the  disturbance  deformation  onto  the  initial  motion  de- 
formation. The  disturbance  deformation  can  be  determined  by  any 
of  the  methods  used  in  elastic  beam  theory. 

a.  Bottasinaafl- ioiaiian 

The  disturbance  deformation  due  to  the  application  of  the 
bending  moment  at  x ° -t  can  be  found  from  Boussinesq’s  solution 
for  an  infinite  beam  until  such  time  as  this  solution  gives  an 
appreciable  disturbance  at  the  fixed  end  of  the  beam.  While  the 
solution  does  not  satisfy  the  fixed  end  boundary  conditions,  it 
yields  an  approximation  to  the  disturbance  deformation  until  suoh 
time  as  the  deflection  and  slope  of  the  beam  at  x = 0 provided 
by  this  solution  become  significant. 

The  Boussinesq  solution  [2]  of  Eq.  (2)  for  a semi  infinite 
beam,  initially  (t  = - oo)  straight  and  at  rest,  and  subject  to 
the  conditions 


Ox^ 


“ 0 , 


= F »(ct) 
8x2  1 


at  x = l , is 


y = 


oo 

% [ Fi(ot 

o 


r 2 2 

5)  cos  ^5*  da 
2a2  2 


(5) 


where  ( = t - x and  F^!(ct) 
ct 

j F^1 (ct)d(ct) • 

-00 


dTct)  ^ence  F^Cot)  ” 
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Now,  since  the  bending  moment  applied  at  the  free  end  of 
the  cantilever  is 


MpCt)  = 


and 


it  is  seen  that 


0 for  t i 0 

9bt  for  t 2.  0 


Mn  = bKT 


d3yr 


ax' 


= o, 


8~yp  _ 


ax' 


2 = Fi*(ct) 


0 for  ct  i 0 

^(ct)  for  ct  2.  0 
c 


at  x = -t. 

It  follows  from  Eq.  (*))  that 
00 


= -2- 


Ac  J 


Cct  - ) cos  ~ do  , 

_ J o * 


H-x)/(2ct)* 


2 a' 


oo 


QyD  = 

ax 


✓nc 


tt-x)/(2ct)* 


(Ct  - lissl  cos  ^ da, 


2ac 


(2 


and 


co 


kd  = 112  = 02 
ax2  Ac 


,2  2 


(6) 


j (ct  “ sin  y*  da.  (7) 
a-x)/(2ct)-? 

From  Eqs.  (3)  and  (7)  it  is  clear  that  the  final  rate  of 
change  of  curvature  of  the  beam  due  to  the  initial  motion  plus 
the  disturbance  motion  is 
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H = if. 5 + 9 |i  - 2S[(-Ux)2Act  3 


(-t-x) 2 Act 


where  8 [(J-x)2Act  ] Is  the  Fresnel  integral, 


\Z2tT 


stiw  dTV 
/n 


Since  S [(-t-x)  ^ Act  3 < .75  for  all  values  of  the  argument,  it 


follows  that 


> 0 


everywhere  along  the  beam  for  all  time.  Thus,  Eq.  (2)  Is  every- 
where valid  and  superposition  of  effects  is  allowable.  Super- 
posing Eqs.  (3)  and  (6),  we  obtain 


y - 2-25t*2  - £ + & 


co 


(•l-x>/(2et)* 


cot  - c-t-iO  V cos  -S?  da. 
2a2  2 


This  Boussinesq  solution  was  calculated  [33  at  t = .01  and 
t = .02  for  the  case  when  t = 3 and  e2  = b/m  = 907.  The  results 
are  shown  in  Table  I.  A plot  of  yD  vs.  x for  t = .02  is  shown 
by  the  dotted  curve  in  Fig,  2. 
b.  joetJiafl. 


This  problem  may  also  bo  solved  by  finding  yD  by  the  well 
known  method  of  normal  modes,  and  then  superposing  this  solution 
onto  the  initial  motion  of  the  beam. 

The  disturbance  deformation,  yD,  is  then  assumed  to  be  of 
the  form 

yD  = ? Xn(x>pn(t)  (8) 

n=l 

where  the  Xn(x)  are  the  normal  modos  for  the  cantilever  beam 
and  the  fn(t)  are  the  normal  coordinates  P+3  of  the  system. 


I 
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Tablo  I 


t = .01 

t = .02 

X 

, 1 
yI 

yD 

7 

1 

yI 

yD 

y 

0 

0 

+.00004 

+.00004 

0 

+.00102 

+.00102 

.15 

-.00024 

-.00005 

-.00029 

+.00026 

+.00173 

+.00199 

.3° 

-.00098 

-.00015 

-.00113 

+.00105 

+.00223 

+.00328 

.45 

-.00219 

-.00019 

-.00238 

+.00236 

+.00245 

+.00481 

.60 

-.00390 

-.00013 

-.00403 

+.00420 

+.00228 

+.00648 

.75 

-.00609 

+.00005 

-.00604 

+.00656 

+.00146 

+.00802 

.90 

-.00878 

+.00029 

-.00849 

+.00945 

+ .00006 

+.00951 

1.0? 

-.01194 

+.00052 

-.01142 

+.01286 

-.00192 

+.01094 

1.20 

-.01560 

+.00062 

-.01498 

+.01680 

-.00438 

+.01242 

1.35 

-.01974 

+.00050 

-.01924 

+.02126 

-.00712 

+.01414 

1.50 

-.02438 

+.00008 

-.02430 

+.02625 

-.00988 

+.OI637 

1.65 

-.02949 

-.00064 

-.03013 

+.03176 

-.01230 

+.01946 

1.80 

-.03510 

-.00156 

-.03666 

+.03780 

-.01399 

+.02381 

1.95 

-.04119 

-.00254 

-.04373 

+.04436 

-.01448 

; +.02988 

2.10 

-.04778 

-.00332 

-.05110 

-.01330 

| +.03015 

2.25 

-.05484 

-.00362 

-.05846 

+.05906 

-.00996 

i +.04910 

2.40 

-.06240 

-.00311 

-.06551 

+.06720 

-.00397 

; +.06323 

2.55 

-.07044 

-.00142 

-.07186 

+.07586 

+.00509 

; +.08095 

2.70 

-.07898 

+.00175 

-.07723 

+.08505 

+.01763 

1 +,10268 

2.85 

-.03799 

+.00668 

-.08131 

+.09476 

+ .03394 

+.12870 

3.00 

-.09750 

+.01355 

-.08395 
■ - — _ . 

+.10500 

+.0^421 

+.15921 

1 
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It  can  easily  toe  shown  toy  the  usual  method  of  finding  normal 
modes*  [5  3 that 


XR(x)  = -j(cos  kjjX  - cosh  knx)(cos  kn-t  + cosh  kn-0  + ( sin  l^x 

- sinh  knx)(sin  kn£  - slnh  kn^J/(m^  sin  knt  slnh  k^ 

(9) 

where  the  k are  obtained  from  the  roovs  of  t he  equation 
n 

cosh  kn-t  cos  kj^  - 1. 

Lagrange's  equations  of  motion  are  simply 

* n + pn2»n  =°n 

in  which  denotes  the  force  corresponding  to  the  coordinate  9^, 

p p if 

and  p * - c^k  • 
n n 

In  order  to  obtain  a force®  .assume  that  a small  increase 

n 

b<Pr  is  given  to  a coordinate  <pn*  The  work  done  toy  the  external 
bending  moment  on  the  end  of  the  beam  is 


V*n  = 9bt6(03T). 


and  so 


. 2k_(sin  kn-t  cosh  k-'t  + sinh  k £ cos  k„£) 

® = 9bt<! n 2 S Q . 

[ -(ml)*  sin  kn£  sinh  kn^  J 


Substituting  ®n  into  the  equations  of  motion  and  taking  into 
consideration  the  homogeneous  initial  conditions  on  yD,  we  find 


* V/hen  the  principal  modes  have  been  normalized  m \ X ‘dx  = 1* 

Jo 
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9n  = 3(t  " pT  Sl“  Pn4) 
Pn  F« 


2kn(sln  kn't  cosh  kn£  + slnh  k^  cos  kn^t) 
sin  kn£  slnh  kn£ 


Finally,  substituting  9n  into  Eq.  (8)  and  superposing  the  initial 
motion  yields  the  solution 


2kn(sln  kn£  cosh  k^  + slnh  k^  co3 
(art)*  sin  kn^t  sinh  kn*t 


>X 


n 


where  the  Xn  are  given  by  Eq.  (9)*  However,  the  convergence  of 
this  expression  is  improved  by  subtracting  out  the  quasi  static 
deflection,  ^ tx2,  due  to  the  application  of  the  bending  moment 
to  the  end  of  the  beam.  Upon  expanding  this  in  terms  of  the 
normal  modes,  it  becomes  clear  that  y may  be  rewritten 

y = 2.2?tx2  - + ^-tx2  - -2k  sin  p t • 

™ 2b  n=l  P 2 n 

n 

{2kn(sin  kn£  cosh  kn£  + sinh  kn£  cos  kn^ 

(mi)'fr  sin  knt  sinh  kn£  Xjl 


This  normal  mode  solution  was  also  calculated  at  t = .01  and 
t - .02  for  the  case  when  t = 3 and  c2  = 907*  Six  terms  of  the 
series  were  used.  The  results  of  these  calculations  are  shown 
in  Table  II,  and  plots  of  yj  vs.  x and  y vs.  x for  t = .02 
are  shown  in  Fig.  3»  A comparison  of  the  disturbance  deformation 
obtained  at  t = .02  by  the  Boussinesq  and  normal  mode  methods  is 
shown  in  Fig.  2, 


0 

..00024 

-.00098 

-.00219 

-.00390 


-.00609  +.00004 
-.00878  +.00031 
-.01194  +.00055 
-.01560  +.00064 
-.01974  +.00049 


.02438 

.02949 

.03510 

.04119 

.04778 

.05484 

.06240 

.07044 

.07898 

.08799 


+.00006 

-.00065 

-.00157 

-.00253 

-.00331 

-.00362 

-.00311 

-.00144 

+.00173 

+.00667 


0 

.00031 

.00116 

.00243 

.00407 

.00605 

.00847 

.01139 

.01496 

.01925 

.02432 

.03014 

.03667 

.04372 

.05109 


+.00026 
+.00105 
+ .00236 
+.00420 

+.00656 

+.00945 

+.01286 

+.01680 

+,02126 

+ .02625 
+.03176 


'D 

0 

+.00029 

+.00091 

+.00149 

+.00172 

+.00135 

+.00027 

-.00154 

-.00399 

-.00685 

-.00980 

-.01242 

-.01424 

-.01474 

-.01346 

-.00997 

-.00385 

+.00^27 

+.01774 

+.03391 


0 

+.00055 

+.00196 

+.00385 

+.00592 


+.00791 

+.00972 

+.01132 

+.01281 

+.01441 

+.01645 

+.OI934 

+.02356 

+.02962 

+.03799 

+.04909 

+.0633^ 

+.08113 

+.10279 

+.12867 


-.097501  +.01357 


+.05397  +.15897  ! 
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3*  Problem  II  --Simply  supported  beam.  Consider  a team  of  uni- 
form cross  section  simply  supported  at  the  ends  x = 0 and  x =£» 
Initially  the  beam  is  at  rest  and  subject  to  an  external  bending 


moment,  Mj  ^ - MQ,  at  both  ends*  At  time  zero  an  external  load 
distribution 


q(x,t)  = qQt  sin 


is  applied  to  the  beam  and  the  resulting  deformation  is  to  be  • 
determined. 


Initially  the  beam  is  entirely  plastic.  If  the  motion  due 
to  the  applied  load  or  disturbance  is  such  that  the  rate  of 

j 

change  of  curvature  of  the  beam  is  everywhere  negative,  for  all 
t,  the  bending  moment-curvature  relationship  of  the  beam  is 


H = - Mq  + bK 

and 

Md  = bKD 

where  MD  and  KD  are  the  additional  bending  moment  and  curvature 
due  to  the  applied  external  load.  Again  the  deformation  of  the 
beam  can  be  determined  by  finding  the  disturbance  deformation  by 
the  normal  modo  method  used  in  elastic  beam  theory,  and  then 
superposing  this  onto  the  initial  deformation. 

It  is  assumed  in  solving  this  problem  that  the  rate  of 
change  of  curvature  of  the  boam  is  everywhere  negative,  for  all 
t.  Once  the  solution  has  beon  obtained,  it  is  roadily  verified 
that  this  assumption  is  valid. 

The  initial  deformation  of  the  boam  is 

yT  = M - x)x 

1 2b 
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and  the  disturbance  deformation  of  the  beam  is  assumed  to  bo  of 


tho  form 


yn  = KMvJ t) 


n Tn 


(10) 


whore  tho  Xfl(x)  are  tho  normal  modes  for  a beam  simply  supported 
at  x = 0 and  x = l • Those  normal  modes  are  readily  found  to  bo 


Xn  = (2/ml r sin  ^ x 


(11) 


Again  wo  havo  as  the  equations  of  motion 


cp  + p m = <X> 

Tn  *n  Tn  n 


(12) 


in  which  <t>n.is  the  force  corresponding  to  the  coordinate  and  1 

Pn2  = • If  a small  increaso  69  n is  given  to  a ooordin- 

ate  9n,  tho  work  dono  by  tho  oxtornal  load  will  bo 

l 


* 

= qQt  sin  ^ (2/mU^sin  DjjL  x6q>n)dx 
« 


and  thus 


<t>  = 0 
n 


for  n/  1 


\ = (2/mt)*tt/2)q0t  . 


Substituting  those  into  tho  equations  of  motion  and  taking  into 
consideration  tho  initial  conditions  of  tho  problom, wo  obtain 


9n  = 0 


for  n / 1 


91  = (2/m£)*tt/2)q( 


(t  - -i-  sin  Pit), 
Pi  * 


Finally,  substituting  Xn  and  9n  into  Eq,  (10)  and  superposing 
tho  initial  doformation,  the  following  solution  for  the  displace* 
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r 

y(x,t)  = ■ -Q(£  - x)x  + — sin  5X  Jt  - JL  sin  p,t 

2b  mp^2  ^ 1 Pi  *■ 


To  ensure  that  this  Is  the  solution  of  the  problem,  It 
remains  only  to  verify  that  ^ < 0 for  0 £ x £ t and  all  t > 0, 
But,  differentiating  y(x,t),  we  find 


1 - COS  Pjt 


which  is  readily  seen  to  bo  negative  for  0 ^ x ^-tand  all  t > 0« 


4.  Problem  III  — A free  boundary  problem.  In  the  two  previous 
problems  the  beam  was  initially  entirely  plastic  and  remained 
entirely  plastic  for  all  time.  No  plastic-rigid  boundaries  were 
present.  The  absence  of  such  boundaries  accounts  in  large  moasure 
for 'the  ease  with  which  these  problems  were  solved.  However,  in 
tho  case  of  a beam  which  is  initially  rigid,  the  application  of 
a load  will  in  general  produce  a moving  plastic-rigid  boundary. 

The  determination  of  this  froo  boundary  and  tho  solution  of  such 
problems  doos  not  appear  simple.  However,  by  tho  inverse  method 
shown  bolow,  it  is  possible  to  find  solutions  to  somo  froo  bound- 
ary problems,  under  the  assumptions  of  Soc.  1.  The  method  con- 
sists essentially  of  determining  what  physical  problems  are 
associated  with  simple  analytic  expressions  satisfying  Eq.  (2). 

Consider  tho  following  solution  of  Eq.  (2): 


y = x 


2t2 


(13) 
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This  expression  satisfies  the  boundary  conditions 

7 * ° and  at  x » C(t)  = (12)*(et)*[K  * ° . 

m * 0 


at  x * 0 < 


That  is,  the  above  solution  of  the  differential  equation  of  motion 
for  plastic  flow  satisfies  fixed  end  conditions  at  x = 0 and  free 
boundary  conditions  at  x - £(t).  As  t increases  from  t = Of  the 
position  of  the  free  boundary  moves  from  the  fixed  end  (x  = 0) 
to  larger  values  of  x.  Moreover,  since  K = 0 for  x *=  £(t),  it 
appears  that  Eq.  (13)  may  be  the  solution  in  the  plastic  region, 
O^xi  £(t),  for  a cantilever  problem  in  which  the  beam  is 
initially  straight.  For  £(t)  < x ^ the  beam  remains  rigid; 
t is  the  length  of  the  beam.  This  solution  will  bo  valid  only 
until  such  time  as  the  plastic-rigid  boundary  reaches  the  free 
end  of  the  beam;  that  is,  for  0 £ t £ l2/c(12)^  , 

The  physical  problem  associated  with  Eq.  (13)  is  now 
assumed  to  be  that  of  a cantilever  beam,  fixed  at  x = 0,  initially 
straight  and  at  rest,  and  subject  at  time  zero  to  a force  F(t) 
and  a bending  momont  M(t),  both  of  which  are  applied  at  the  free 
end,  x -l  • 

F(t)  and  M(t)  are  determined  from  the  equations  of  motion 
for  the  rigid  portion  of  the  beam,  i.e., 


F(t)  + b(y,M)x=^  =j 


(R 


l 

* 

ndR  + ymdR 

« 

K 


1 1 

f 

M(t)  + F(t)  (*t  - O = M#  + 

w 


(R 


O 


j(R 


(R  -OymdR. 
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It  follows,  then,  that 

F(t)  = - b(y"*)  _ + - O(y)  r + m lA."  -0^(9’ 

x=<  x=£  2 1 

M(t)  - - FCtH't  - O + M0  + m Ii_=JLl3(e)v_r  + a JLLl. 

° 3 x- ^ 2 


’•/here 


5(t)  = (12)*(ct)* 


= 6(12)*ct 


(e)x=C  = ( -2|JC-  + 25  ♦ «>2  ^L.  = 6(l2)±(ct)* 

^ dt2ax  ax2at  ax3  x ^ 


ax2at 


(y"‘). 


- . 


For  any  given  set  of  parameters  a check  must  be  made  to  make 
sure  that  M(x)  < M0  for  5 < x £ Z • That  is,  it  must  be  verified 
that  the  yield  limit  is  not  exceeded  in  that  portion  of  the  beam 
which  is  assumed  to  be  rigid.  M(x)  is  readily  obtained  from  the 
follov/ing  equilibrium  equation: 


M(x)  = M(t)  + F(t)(J  -x)  -I  (R  - x) j(R  - K )©|mdR 


my(R  - x)dR  • 


As  a specific  example,  the  following  practical  values  of 
the  parameters  were  chosen:  Z = 5 ft.,  M0  = 17  ft.-lbs.,  m = 
1.2  lbs. /ft.,  c2  = 907  ft.^/sec.2  • H(x)  was  calculated  for 
5 £ Z at  t = 0,  .005,  ;02,  .05,  .10,  .20  secs.  From  the 
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results  of  these  calculations,  v/hich  are  shown  in  Fig*  4,  it 
appears  that  for  this  set  of  parameters  M(x)  < M0  for  K < x £ l. 
It  does  not  seem  feasible  to  make  a rigorous  analytic  check  of 
this  fact* 

5*  Conclusion.  The  examples  troated  in  this  report  serve  to 
illustrate  the  fact  that,  under  certain  conditions,  tho  plastic 
doformation  of  beams  subject  to  dynamic  transverse  loading  can 
be  determined  by  the  same  techniques  used  in  solving  elastic 
beam  problems*  This  is  due  to  the  fact  that,  under  thG  assump- 
tions stated  in  Sec.  1,  tho  equation  of  motion  for  a beam  under- 
going largo  plastic  deformation  is  of  tho  same  linear  form  as 
the  equation  of  motion  for  an  elastic  beam* 

Problems  I and  II  wero  readily  solved,  largoly  due  to  tho 
fact  that  in  oach  case  the  beam  was  initially  entirely  plastic 
and  remained  entirely  plcstic  for  all  timo.  If  the  beam  is 
initially  rigid,  an  applied  load  will,  in  general,  produce  a 
moving  plastic-rigid  boundary*  Tho  presenco  of  such  froe  bound- 
aries makes  those  problems  more  difficult  to  handle.  The  same 
difficulty  arises,  of  course,  in  solving  probloms  involving  the 
unloading  of  external  forces,  rathor  than  tho  loading  typo  of 
problems  troated  hero. 

As  shown  by  Problem  III,  it  is  possible  to  find  solutions 
to  some  froo  boundary  problems  under  tho  assumptions  of  Sec*  1* 
Tho  method  of  solution  is  ossontially  an  inverse  procoss  wheroby 
a simplo  analytic  cxprosMon  satisfying  the  differential  equation 
of  motion  for  plastic  flow  is  considered,  and  the  physical  prob- 
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lem  associated  with  this  solution  is  then  determined.  The  general 
problem  of  a beam  subject  to  prescribed  external  forces  giving 
rise  to  free  boundaries  is  not  so  readily  solved. 

It  was  hoped  that  a numerical  technique  based  on  Hartroe*s 
mothod  [6]  might  be  devised  for  locating  these  free  plastic-rigid 
boundaries  at  any  given  time.  While  there  was  no  free  boundary 
to  be  determined,  Hartree  rs  method  was  first  tried  on  Problem  I 
of  this  report.  However,  it  was  found  that,  due  to  the  prosonce 
of  the  fourth  derivative  in  tho  equation  of  motion  of  tho  beam, 
the  results  of  a finite  difference  technique  to  solve  the  ordinary 
differential  equation  resulting  from  Hartroofs  mothod  were  poor* 

In  the  caso  of  a free  boundary  problem,  tho  satisfaction  of  tho 
boundary  condition,  ^ = 0,  would  introduce  further  high  deriva- 
tives, and  it  was  virtually  impossible  to  retain  sufficient 
accuracy  using  finito  difference  to  complete  such  a problem* 
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